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Chapter 1 
INTRODUCTION 

Although the elastic finite-displacement response of shallow cylin- 
drical ly curved panels has been subjected to extensive consideration in 
the past fifty years, very little work has been done on the behavior of 
such panels fabricated from orthotropic materials. The utilization of 
composite materials in primary structural components in aircraft appli- 
cations is steadily increasing. The current trend is the implementation 
of composite materials into the fuselage structure. Figure 1 shows how 
the fuselage skin may be represented by a shallow cyl i ndri cal ly curved 
panel . 

One concern that must be addressed involves the crashworthiness of 
a composite fuselage. Crashworthiness studies imply high speed, large 
deflection behavior of structural components. It is the intent of this 
study to lay the groundwork for understanding the dynamic, elastic, 
finite-displacement behavior of shallow cyl indrical ly curved panels 
fabricated from a fibrous composite material. This will be accomplished 
by presenting a complete examination, both theoretical and experimental, 
of the static behavior of such a structure subjected to a line load 
acting along a cylinder generator. The line load acts toward the center 
of curvature of the panel. This examination of static behavior is a 
necessary step before any attempts at determining dynamic behavior are 
made. 
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FIG. 1 Model of fuselage skin 
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Previous Work 

In 1929, Biezeno [1] investigated the response of a shallow circu- 
lar arch with pinned ends at a fixed distance apart subjected to a 
concentrated load acting at midspan. Using two solution methods, one a 
shallow curved-beam solution and the other a series solution, Biezeno 
presented a detailed analysis of the symmetric response. Herein, symme- 
try is defined with respect to displacement and the midspan location on 
the curved member. Figure 2a illustrates what is referred to as a 
symmetric response while Fig. 2b ilVustrates an asymmetric response. No 
results were presented by Biezeno concerning the asymmetric response 
except to mention its existence and its importance to loss of stability 
of the symmetric equilibrium configuration. Following the approach of 
Biezeno, Fung and Kaplan [2] conducted a complete investigation into the 
behavior of shallow arches, specifically the case of an initially sinu- 
soidal arch, using a series solution. They investigated the asymmetric 
displacement mode in greater detail than did Biezeno. It was found that 
on the primary load-deflection equilibrium path, corresponding to sym- 
metric deformations, asymmetric equilibrium states occurred at unstable 
bifurcation points. Schreyer and Masur [3] performed a similar analysis 
on the circular arch with clamped ends. In addition to determining the 
equilibrium behavior of the structure, they also performed an in-depth 
analysis of the stability of these equilibrium states. In their study 
the initial rise-to-span ratio determined whether the buckling of the 
arch occurred at a limit point on the primary path or whether it 
occurred at the bifurcation point on the primary path. These first 
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(a) Symmetric 
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Asymmetric 


FIG, 2 Symmetric and asymmetric response 
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studies into the finite-deformation behavior of the shallow arch were 
restricted to isotropic and homogeneous materials with linear elastic 
behavior and did not consider load eccentricity or geometric imperfec- 
tions. 

In recent years it became obvious that a theory which did not take 
into account initial imperfections of the structure would not be able to 
accurately model a real structure. Due to the imperfection sensitivity 
of some arch structures, it was at times difficult to experimentally 
verify the analyses put forward. The buckling loads, as well as the 
entire character of the response, are affected by the imperfections 
present. Schreyer [4] presented an analysis of the effects of initial 
imperfections on the response of shallow circular arches with clamped 
ends. In the study, imperfections in the geometry of the arch, which 
had the same form as the buckling mode, and eccentricities in the load 
position were both examined. Both types of imperfections were seen to 
cause a decrease in the buckling load for certain ranges in geometry. 
This implied that the arch is an imperfection-sensitive structure. 

Plaut [5] performed a more thorough investigation into the effects of 
load position on the magnitude of the buckling load. Examining the 
pinned-ended shallow circular arch, a number of different behaviors were 
observed due to load offsets. However, only those arches which exhibi- 
ted loss of stability at bifurcation points when the load was exactly at 
midspan were seen to be sensitive to load offset. In such an arch a 
small perturbation in the load position could cause significant reduc- 
tions in the predicted buckling loads. Cheung and Babcock [6] performed 
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an experimental investigation into the effects of load offset on the 
static buckling load of shallow circular arches. Using aluminum arches 
with clamped ends, buckling loads were measured for load offsets of up 
to 6% of the arc length. The results of the experiment show significant 
reductions in the magnitude of the buckling load with increasing load 
offset. Taking into account such load-position effects considerably 
improved the understanding of the response of the shallow arch. How- 
ever, the vast majority of the work has concentrated on metallic 
arches. With composite materials playing an increasingly important role 
in the aerospace industry, the effects of orthotropic materials on the 
response of the shallow arch was still needed. 

Some work on the behavior of orthotropic shallow arches has been 
done. Nash and Hsu [7] analyzed the structural behavior of a shallow 
circular arch fabricated from a composite material with the fiber direc- 
tion normal to the middle surface. In their theory, Nash and Hsu inclu- 
ded the effects of transverse shear stresses and transverse normal 
stresses. The qualitative behavior of their orthotropic arch was essen- 
tially the same as for an isotropic arch. As is the case with an iso- 
tropic arch, they found the behavior of the arch was controlled by a 
geometric parameter. In addition to the analytic considerations, they 
also conducted experiments into the load-displacement relationship up to 
the buckling load. In the results obtained, a fair agreement between 
theory and experiment was seen. Marshall, Rhodes, and Banks [8] conduc- 
ted experiments into the behavior of shallow curved composite panels, 
also finding good agreement between theory and experiment. These exper- 
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iments were performed on spherically curved glass-epoxy panels under 
load control. 

Objective of Present Work 

The study discussed in this report is intended to extend the pres- 
ent knowledge of the static, large deflection, behavior of orthotropic 
shallow curved panels. To accomplish this, an exact solution is presen- 
ted for a simply supported cylindrical panel with a circular cross 
section and subjected to a line loading along a cylinder generator. The 
orthotropic axes are parallel to the generator and circumferential 
directions. Derivation of the equations is presented in Chapter 2. The 
governing equations are derived using laminated plate theory, nonlinear 
strain-displacement relations, and applying variational principles. 
Eccentricities, or offsets, in load position relative to midspan are 
taken into account. Due to the geometry of the problem under considera- 
tion, only a one-dimensional analysis is developed. The one independent 
variable is related to the ci rcumferential arc length along the panel. 
Because of the one-dimensional nature of the problem, ordinary differen- 
tial equations govern the behavior of the panel. These equations are 
solved in closed form. The aim of this work is to investigate the 
effects of orthotropic material properties on the behavior of a 
shallow curved panel. In past work, as has been discussed, the behav- 
ior of the structure is assumed to be controlled solely by the geome- 
try. For an orthotropic curved panel, the material properties may also 
affect the response. This idea is investigated here. In addition, a 
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load-position sensitivity analysis is presented along with predictions 
of load vs. displacement, load vs. thrust, and load vs. surface 
strains. These predictions are presented in Chapter 3. In order to 
verify the analysis, experimental results of displacement-controlled 
tests performed on graphite-epoxy curved panels are presented. Chapter 
4 describes the experimental set-up while Chapter 5 presents experimen- 
tal results. Since the tests are run under displacement control, a 
greater portion of the equilibrium path is stable and predictions of the 
buckling load under load control may be made with greater confidence. 
Also, experimental observation of the post-buckling mode, unstable under 
load control, may be made using displacement control. By comparing 
analytic and experimental results, the effectiveness of laminated plate 
theory, in combination with an analysis accounting for large elastic 
deformations, in predicting the structural response will be shown. In 
addition, the influence of panel geometry and panel material properties 
on the response will be illustrated. Finally, Chapter 6 presents con- 
clusions and recommendations for future work. 


Chapter 2 

DERIVATION OF GOVERNING EQUATIONS 

Nomenclature 

The coordinate system and notation used for defining the cylindric- 
ally curved panel is given in fig. 3. The panel has a radius of curva- 
ture R. The radius is measured from the center of curvature, denoted by 
0, to the geometric middle surface. The semi -opening angle is denoted 
by 0. The independent spatial variables are represented by 9 and x. 

The angle 0 is measured positive clockwise from the midspan vertical and 
the generator coordinate x is measured positive from front to back on 
the panel. In addition x = 0 at the front. Position through the thick- 
ness is given by z, with z being measured positive away from the center 
of curvature. Furthermore, z = 0 at the middle surface. The total panel 
thickness is denoted by h. The panel is assumed to be simply supported 
with a fixed distance L $ between supports. The length of the panel in 
the x direction is L x . The panel rise, H, is the vertical distance, 
measured at midspan, from the line of supports to the middle surface. 

The radial displacement of the middle surface is given by w. The 
radial displacement is positive toward the center of curvature. The 
tangential displacement, represented by u, is positive toward the 
right. The line load intensity, assumed to be independent of the x 
coordinate, is denoted by P and is assumed to act on the middle surface 
at an arbitrary spatial position denoted by "0. 

Due to the one-diinensional character of the geometry and loading, 
the displacements are assumed to be functions only of the spatial vari- 
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FIG. 3 Cylinder geometry and nomenclature 
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able 9. Thus in the analysis, w = w(0) and u = u(0). Other assumptions 
to be used in the analysis include the assumptions of a linear elastic 
material behavior and linear strain variation through the thickness. 
Also, the panel is assumed to be thin with the transverse normal stress 
negligible in comparison with the inplane stresses. Transverse shear 
deformations are neglected. 

Kinematics of Deformation 

The kinematics of the middle surface deformation are shown in fig. 
4. The applied load causes each point of the panel middle surface to 
displace by an amount t and the corresponding cross section to rotate by 
an amount Q. The unit vectors to be used for the undeformed geometry 

a 

are i p , normal to the middle surface and directed toward the center of 
curvature, and the tangential unit vector i^. For the deformed config- 

A A 

uration the radial and tangential unit vectors are e r and e^, respect- 
ively. The position vector of a point on the undeformed middle surface 
is given by, 

* = - Ri r - (1) 

The position vector for the same point after deformation is denoted by, 

r = It + t. (2) 

The displacement vector is 

t = w(S)i p + u (S ) i t , (3) 

where the undeformed arc-length is 

S = R0. (4) 

The length of a differential element of the undeformed arc length 
is denoted by dS. The length of this differential element of arc length 



FIG. 4 Kinematics of middle surface deformations 
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in the deformed configuration is represented by dS*. By definition the 
extensional strain of the middle surface is. 


o _ dS*-dS 
e 9 " dS * 


( 5 ) 


To derive the governing equations in terms of displacements, it is 
necessary to find the relation between the arc lengths and displace- 
ments, This relation can be determined by realizing, 


dS* = 


| dS l dS * 


( 6 ) 


From eq. 2 the change of the deformed position vector with respect to 
the undeformed arc length is. 


dr dft A di 

<15 ^ dS’ 


( 7 ) 


where. 



(8) 


and. 



( 9 ) 


To obtain eqs. 8 and 9 use was made of the relations 
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di r It di t It 

35 R \ and dS - " + R V* ( 10 )’ 

Realizing that + represents the initial curvature, < 0 , the displace- 
ment gradients, r and r. , are defined as, 
n l 

r n = 7K + K o u and r t = 3? " K o w ( 12 )» ( 13 ) 

so that eq, 7 becomes. 


dr 

33 


= r nV 


+ (1 + r t )i t . 


(14) 


From eqs. 6 and 14 the deformed arc length is related to the undeformed 
arc length by the relation, 

dS* = [(1 + r t ) 2 + r 2 ] 1/2 ds. (15) 

The extensional strain, e°, is now calculated to be, 

u 

e 0 = + r t )2 + r nJ 1/2 - 1 • 06) 

Squaring the first term in the radical and applying the binomial expan- 
sion leads to. 


c e “ r t + ' /2 r n + °( r3 >- 


(17) 
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By assuming the displacement gradients remain small, i.e.. 


0 < I r t I « 1 and 0 < |rj « 1, (18), (19) 


terms of third order and higher are dropped from the strain expres- 
sion. This results in the approximate strain-displacement relationship, 


o _ du 
e 0 dS 


K 0 W 


1 


'HS + 


2 'dS 


v> 


( 20 ) 


The configuration of the deformed middle surface is specified by 
the position vector r(S*). Let k(S*) represent the curvature of the 
normal section of this surface perpendicular to the x axis. Then the 
differential geometry of this normal section requires 


dr 

dS* 

A 

= e t 

(21) 

A 

de 

A 


ii 

-i 

- K e t 

(22) 

A 

de t 

A 


d3* 

= <e p 

(23) 


These expressions are analogous to eqs. 8, 10, and 11 for the undeformed 
normal section of the middle surface. The unit tangent vector, eq. 21, 
is related to the displacement gradients, eqs. 12 and 13, by using the 
chain rule to relate S* to S, and then using eqs. 5 and 14. The result 
is 
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e t = [ r nV + ( 1+r t) i t^ /(1+E e ) * ( 24 ) 

Now define the clockwise angle of rotation of the normal section due to 
deformation as Q; see fig. 4. Thus 


A /V 

= si n (£2 ) i r + cos(Q)i t , 

(25) 

A A 

= cos (S2 ) i p - si n (£2)i t . 

(26) 


Comparing eqs. 24 and 25, the trignometric functions of the rotation 
angle are 


sin (£2) 

n 

(27) 


1+e 9 



1+r t 


cos ( q ) 

u 

r» 

(28) 


1+e 0 


tan(Q) 

r n 
' 1+IV 

(29) 


Consider the deformation of a surface parallel to the middle sur- 
face specified by a nonzero constant value of z. Let dS z and dS z * 
designate the undeformed and deformed differential arc-lengths, respect- 
ively, of the normal section to the parallel surface. Consistent with 
the definition of the middle surface strain, eq. 5, the circumferential 
strain of the parallel surface is defined as. 
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( 30 ) 


A material point located by the position vector $ z on the undeformed 
parallel surface is located by position vector r on the deformed par- 
allel surface. It is assumed that normals to the undeformed middle 
surface remain normal and undeformed with respect to the middle surface 
during deformation (Love-Ki rchhoff hypothesis). Thus these position 
vectors are 


Using eqs. 8, 10, 21, and 22, the differentials of these vectors with z 
held constant are 


= - (R+z)i r . 


. -* 

and r z = r - ze r . 


(31), (32) 


dlT z = - (l+ZK Q )dS i t , and dr z = - (l+z<)dS* e^. (33), (34) 


The magnitudes of these differential vectors are 


dS z = (l+z< 0 )dS, and dS* = (l+ZK)dS* . (35), (36) 


Substitute these into eq. 30, and use eq. 5, to get 



( 37 ) 


For thin shells. 
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0 < |zk 0 | < |h< Q | < < | . 


( 38 ) 


Assume the shell is thin, and expand the denominator of eq. 37 in a 
power series in z using the binomial series to get 


Consequently, for thin shells and small strains the parallel surface 
strain is approximated by. 


This result shows the parallel surface strain is determined by the 
middle surface strain and its change in curvature. 

To complete the kinematics, the change in curvature of the middle 
surface is needed. According to eq. 23 the curvature of the normal 
section of the deformed middle surface is 


e 


e - 4 + (i+4)[z(k-< 0 ) + o(z 2 )] 


(39) 



( 40 ) 


K = V dS* 


( 41 ) 


Using the chain rule and eq. 5 this becomes 



( 42 ) 
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From eqs. 10, 11, 25, and 26, the derivative of the deformed tangent 
unit vector with respect to the undeformed arc length is. 


A 

dG t f 
35" = ( K o 



(43) 


Substitute this into eq. 41 to get 




(44) 


<>S> 


Using eqs. 27 and 28, the derivative of the rotation with respect to the 
undeformed arc length is. 


dQ 

35 


1 dr n dr t 

~ o tcos Q 33“ " sin 33“ ] 
+e e 


(45) 


For small strains and rotations this is approximated by 


dfi _ dr n 
dS dS * 


(46) 


The implication of the approximation in eq. 46 is important. If the 
rotations are assumed to be small, and if the load-deformation behavior 
of a curved panel from no-load to complete snap-through is to be invest- 
igated, then the initial curvature of the panel must be small. This 
implies the work here is limited to shallow shell theory. The curva- 
ture, eq. 44, can now be expressed as 


K 


(47) 
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dr n 

de- 


noting that the small strain assumption leads to the approximation. 


1 


+ £, 


1. 


(48) 


Substituting for the displacement gradient, the curvature-displacement 
relationship can be determined from eq, 47 as 


K = 




(49) 


From shallow shell theory the tangential displacement u is assumed 

dw 

to be small relative to ^ so that the normal displacement gradient 
becomes , 



(50) 


Thus for shallow panels, the strain-displacement and curvature-displace- 
ment relations are, 


(51) 

(52) 


The change in middle surface curvature, K, is defined to be 

K = k - k q = d^w/dS^ 


du 

dS 


< Q w 


1 ,dw N 2 


d^w 

7 ? 


( 53 ) 
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Total Potential Energy 

As previously stated, the governing equations will be derived using 
variational principles. To do so it is necessary to form the total 
potential energy of the system. By definition the total potential 
energy is, 

* = U + V (54) 


where U is the strain energy stored in the structure and V is the poten- 
tial of the loading. The strain energy density, U^, is defined as. 


U 


d 


1 

2 


a. . e. . 
1J 1J 


(55) 


Since transverse shear strains and the transverse normal stress are 
assumed to be negligible in comparison to the inplane stresses and 
strains, the strain energy density reduces to. 


u d ■ 1 K E x + "xe'xe * Vet 


(56) 


The previously mentioned one-dimensional nature of the geometry suggests 
that, 

^ = 0 and u x = 0, (57) 


where u x is the displacement in the x direction. Using these assump- 
tions the strain energy density is reduced further to. 
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U d = 2 cr Q e 0* (58) 

Since the variables are independent of x, the total strain energy, using 
eq. 58, is. 


L x 

U = y- f f a 0 e e dSdz , (59) 

where L x is the length of the panel in the x direction. 

The constitutive law for an orthotropic material in the state of 
plane stress described here is. 



where the TJ matrix is the transformed reduced stiffness matrix. Equa- 
tion 60 is written in structural coordinates. For a more complete 
treatment of orthotropic constitutive relations see [9]. The stress in 
the 9 di rection is. 


°9 ^22 e 9* 


( 61 ) 


Substituting this relation into eq. 59 yields, 
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u = jr I J Q 22 e e dSdz * 

Strain away from the middle surface, from eq, 40, is written i 
the middle surface strain and curvature change as, 

o 

e e - e 0 + zK . 

With this expression the constitutive law is, 

a 0 = ^ 22^9 + zK ^* 

Integrating this expression with respect to z yields, 


N e = A 22 e 9 + B 22 K * 


the force resultant, N 0 , being defined as. 


N 0 = / a 0 dz . 
z 


The extensional stiffness, A 22 , and bending-stretching coupling 
ness, B 22 , are defined by. 


A 22 = ^22^ z k " Z k-1 ^ 

B 22 = 7 |( J 1 ^22^k ( z k ‘ z k-l ) ’ 


(62) 

terms of 

(63) 

(64) 

(65) 

( 66 ) 
stiff- 

(67) 




(68) 
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where k is the lamina number and and (z^ > z^) define the 

location of the kth lamina in the laminate. Multiplying the constitu- 
tive law of eq. 64 by z and integrating over z yield the expression for 
the moment resultant, M Q , 


where 


M 0 = B 22 e e + D 22 K ’ 


M q = a Q zdz. 


(69) 


(70) 


The bending stiffness, is defined as. 


°22 ' 1 |< f 1 ^22^k^ z k " z k-l ^ 


(71) 


Since only symmetric lay-ups are to be examined, the stiffness is 
set equal to zero in the analysis. 

Substituting the strain expression of eq. 63 into eq. 62 and inte- 
grating over z yields. 


U = f {^22^ e 9^ + ^22^^^ * (72) 

* 

Taking the first variation of the strain energy and substituting the 
constitutive relations of eqs. 65 and 69 into eq. 72 leads to, 


«U * L I 

S 


fN 0 6e° + M 0 6K}dS . 


(73) 


25 


From eqs. 51 and 53 the variations of strain and curvature change are 
written in terms of the displacement functions as. 


o ,du. . dw ,dw. 

6£ 0 6 ^dS^ " k o 6w + dS “^dS^ 


6K 


(^ 5 ) 

dS^ 


(74) 

(75) 


Upon substitution of these relationships into eq. 73 and adding the 
variation of the load potential, the first variation of the total poten- 
tial energy is, 


6n = 


L x Jg t N e 6 ^ " Vo 6w + N e a? 


2 

+ M 0 6(^)}dS - PL x 6w(S) 
dS 


(76) 


The arc-wise location where the line load acts, "S, is given by 


= Re . 


(77) 


Furthermore, one-half the arc length of the undeformed panel, S 0 , is 
given by 

S 0 - RP . 


(78) 
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The theorem of stationary total potential energy states that the 
first variation of the potential energy is zero for equilibrium config- 
urations. By setting 6 ti to zero in eq. 76 and integrating the first, 
third and fourth terms by parts, with the integration taking place over 
the intervals -S 0 < S < S’ and S> < S < S 0 , the following relationship 
results : 


s s n dN 0 . .2 .2 m dN. 

f_ s l{- k o N 0 - d$ dS - N 0 + ds 2 l 6w - dS 6u ^ dS 
o 


* N 9 IS H? S _ * N 9 dS 


,dw» il> 


,dW,, S 0 dM 9 


dw 6w| 0 
■? 

S. 


dM, 


+ M e 6 ©l-S + M 9 6 ©U U " cR 2- 6w l-S " ’S’T *1” < 79 ) 

O S 0 S 


- P6w(3‘) = 0 . 


It was necessary to break the integration into the intervals due to the 
discontinuous nature of the loading. Since the variations on the dis- 
placement are assumed to be arbitrary, the following governing equations 
and boundary conditions are obtained: 


dN 0 

3S“ 


k O N 9 - 


^idw N 

dS dS ‘ "9 


d 2 w 

dS 2 


d 2 M, 


dS 


= 0 


(80) 

(81) 


over the intervals - S Q < S < T and 3T < S < S Q . At S = -S 0 and S = S 0 ; 
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either N = 0 or u is known, (82) 

0 

Hu, dM fl 

either N 0 = 0 or w is known, (83) 

either M Q = 0 or is known. (84) 

The transition conditions at S = TT , 


u is continuous. 

(85) 

w is continuous. 

(86) 

dw 

is continuous, 

(87) 

M 0 is continuous, 

(88) 

N is continuous, 

(89) 

(5*) dM fT) 
s - ds - p = 0 - 

(90) 


Equation 90 indicates that the shear force experiences a jump at S = S. 

An additional equation can be obtained by integrating both sides of 
eq, 51 over the arc length S Q < S < S Q . Using eq. 65, recalling B 22 ~ 

0, eq. 51 assumes the form 


N e 

^22 


du 

S5 


< 0 w + 


1 (dw-i2 

1 W ‘ 


(91) 
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From eq. 80, N. is constant with S. Therefore integration of eq. 91 
results in 


"e " VT A K“ * 7 (^ 2 ) dS • 
0 - s o 

(92) 

where u vanishes at the fixed ends. Using eq. 80 in eq. 

81, results in 

„ „ d 2 w d2 "e „ 

- “h". - "e ;£Z + 

(93) 

Equations 92 and 93 are the equations governing the 

deformation of 

the panel. For the panel configuration in question (see 

fig. 3) the 

boundary conditions are. 


u = 0 

(94) 

w = 0 

(95) 

M 0= 0 

(96) 

at S = + S Q . Transition conditions, applied at S = S, are. 

u (S ) = ufS 4 ) 

(97) 

w(T) = wCT 1 " ) 

(98) 

afon • 

(99) 

M q (T) = M^) 

(100) 

N e (T) = N Q (T f ) 

(101) 

dM 0 (r) dM 0 (T f ) n 
dS " dS - P - 0 . 

(102) 
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In order to facilitate the analysis, the governing equations are 
nondimensional ized. Since 

dS = Rde , (103) 

then 


1 d' 


.~n n n , n * 

dS R d0 


(104) 


The nondimensional quantities are defined as follows 

N 


T) = 


4 > = 


’e 

2 

(105) 

a 22 p 


w 

o 

(106) 

Rp^ 


u 

73 

(107) 

Rp 


0 

" P 

(108) 

P 

. „3 

(109) 
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Using this nondimensional scheme along with the constitutive relations 
of eqs. 65 and 69, the governing equations can be rewritten as. 


44-u4*i).o 


~T — 4 

dr 


dr 


( 110 ) 


and 


1 


11 " " \ ^ " J * 


( 111 ) 
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where X is a nondirnensi onal parameter defined as. 


X 


4 


A 22 

D 



22 


( 112 ) 


Equation 110 corresponds to eq. 93 while eq. Ill corresponds to eq, 92. 

By using the definitions of the bending and extensional stiffnesses 
for an isotropic material, the parameter X is, 

\ 4 = 48(£) 2 . (113) 


It is seen that the material dependence is eliminated and X is a func- 
tion of geometry only. In contrast, with an orthotropic material X is 
also a function of the material properties. It should be emphasized 
that the solution to eqs. 110 and 111 depend only on the parameter X. 
Since, for an orthotropic material, X incorporates both material and 
geometric properties of the panel, all panel characteristics influence 
this single parameter. 

The nondimensional boundary conditions and transition conditions, 
using eqs. 105-109, are, 

u « 0 at r * i 1 (114) 

<(> = 0atr=±l (115) 

= 0 at r = i 1 
dr 


( 116 ) 
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u(T) = u(f f ) 

(117) 

♦(T) = ♦ (f f ) 

(118) 


(119) 

d 2 <f(r") _ d 2 ^) 
dr 2 dr 2 

(120) 

n(T) = n(f + ) 

(121) 

d 3 <t>(r~). = d\(p. + pX *. 

(122) 


dr dr 


Chapter 3 


SOLUTION TO GOVERNING EQUATION AND NUMERICAL RESULTS 


The governing differential equation developed in the previous 
section, eq. 110, appears to be a linear ordinary differential equation 
with T| being an unknown constant. This is not the case since t) is 
actually a nonlinear function of the displacement. However, linear 
analysis may be used to find the functional form of the radial 
displacement, <t>(r). The nonlinearity of the problem manifests itself 
through the second governing equation, eq. 111. Since the functional 
form of <t>(r) is established by eq. 110, eq. Ill results in a transcen- 
dental equation for n. Since only compressive middle surface strains 
(i.e., N 0 < 0) are of interest, the thrust parameter q is defined such 
that. 


Substituting this expression into the governing equations yields. 


Equation 124 can be solved in closed form. However, the form of the 
solution depends on the value of q. The solutions for various values of 
q are discussed below. 



(123) 



(124) 


(125) 
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Midspan loading: q / 0 and q l -jp, n = 1, 3, 5, ... 

With the load acting radially at midspan, T = 0,and with q £ 0 and 
q / -jp, wl1ere n = 1, 3, 5, ...» the general solution to eq. 124 for the 
radial displacement is 

4>(r) = A x + A 2 r - r 2 + A 3 cos(qr) + A 4 sin(qr) (126) 

in each of the intervals -1 < r < 0 and 0 < r < 1. Applying the bound- 
ary conditions, eqs. 114-116, on each interval, and the transition 
conditions, eqs. 117-122, at r = 0 yields, 

4>(r) = b Sin(qr) + - 2 (1 - ^osfqj ) + \ U “ ^ ) 

4 q (127) 
- {1 - |r| tan(q)cos(qr) + j s1n(q|r| )} . 

2q H H 

This solution is valid over both intervals. In this expression b repre- 
sents the magnitude of the asymmetric displacement response. The fol- 
lowing conditions on b are determined from the boundary condition 
requiring the moment to vanish at the edges: 

if q / nit then b = 0 
if q = nn then b / 0, 

where n = 1, 2, 3 To complete the solution process the func- 

tional form of 4) (r ) , eq. 127, is substituted into eq. 125. This relates 
the thrust parameter q to the nondimensional applied load, p. After 
some algebra this relationship is written as, 
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A£ 2 + BC + C = 0 , 


(128) 


where A, B, and C are functions of q and X, and 



(129) 


Appendix A contains the functional form of A, B, and C for all cases 
considered here. For a given panel, X is specified. In eq. 128, by 
specifying the thrust parameter q, a quadratic equation can be solved to 
obtain the corresponding load, p. This is in opposition to solving eq. 
128 as a transcendental equation for the thrust parameter if loads are 
first specified. In this solution procedure it is important to recog- 
nize that eq. 128 may have two real roots, one real root, or no real 
roots, depending on the values of X and q. Each real root corresponds 
to a separate equilibrium configuration. Having found the relationship 
between load and thrust, the load-displacement relationship follows from 
eq. 127. 

For q = nn, b / 0, there is an equation similar to eq. 128 which 
relates the applied load to the magnitude of the asymmetric deformation, 
b. It also originates from eq. 125. This relationship is given by the 
expression, 


Here A, B, and D are constants depending only on the value of n. Con- 
stant C depends on n and X. Again these constants are expressed in 
Appendix A. Real solutions for the asymmetric amplitude b require the 
right-hand-side of eq. 130 to be nonnegative. Consequently, not all 


Ob 2 = AC 2 + BC + C. 


(130) 
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integer values of n are permissible for a specific value of Notice 
that two asymmetric configurations are determined from eq. 130 at given 
load if the right-hand-side is positive. 

Midspan Loading: q = 0 

When the thrust is equal to zero eq. 127 is not the correct func- 
tional form of 4> (r) . Setting q = 0 in the governing equation, eq. 110, 
a new functional form of 4>(r) is obtained. Applying the boundary condi- 
tions (eqs. 114-116) and the transition conditions (eqs. 117-122) the 
functional form of the radial displacement is, 

*(r) = 5 + * (131) 

By substituting eq. 131 into eq. 125 the relationship, 

l Z - ^ l = 0 , (132) 

is obtained for determining the loads at which the thrust parameter, and 

thus the middle surface thrust, N a , vanishes. This occurs at either 
25 

£ = 0 or £ = — ^ . The £ = 0 solution states, as expected, that there 
is no thrust in the panel if there is no load. However, the nonzero 
solution to eq. 132 indicates that there is a loaded and deformed state 
of the panel for which there is no thrust in the panel. In this special 
situation the load is resisted by bending only. 
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Midspan Loading: q = , n = 1, 3, 5, ... 

The functional form for <f> (r) in eq. 127 is incorrect when q = -jp, 
n = 1, 3, 5, ... . In this case the linear inhomogeneous mathematical 
problem for <t> (r ) , given by eq. 124, the boundary conditions, and the 
transition conditions, does not have a unique solution. The mathemati- 
cal reason for this is found in Fredholm's alternative theorem [10, pg. 
547]. The associated adjoint homogeneous problem for <t> (r) has nontriv- 
ial solutions (eigenfunctions) associated with these values of q (eigen- 
values). Consequently the conditions for a solution to 4> (r ) require, 

+ 1 n 

l 4> m (0) = / q‘ <t> m (r)dr , (133) 

where d> (r) = A cosfS^-), m = 1, 3, 5, ... are the eigenfunctions of 
m ' ' m v 2 1 

the associated homogeneous problem. (Equation 133 is equivalent to the 
"inhomogeneity" being orthogonal to the eigenfunctions.) Equation 133 
determines the load for the given values of the thrust. The displace- 
ment functions which result from this are, 

, o * 2sin(^) 

4> n (r) = \ (l - r) + A n cos(qr) + ( l r | - l) 

nit 

(134) 

* (|/ U - Sln(j£) s1n(q|r| )} , 


n = 1, 3, 5, 


n max • 


* 

The values of A n 


are found by substituting eq. 


134 into eq. 125. This yields 
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A(A*) 2 + BA* + C = 0 , (135) 

where A, B, and C are constants depending on the values of n and x. 
Appendix A presents the coefficients A, B, C. The solution for <t> is not 
unique in this case. The finite number of these solutions, n max , 
follows from the requirement of real solutions to eq. 135. 

Eccentric Loading: q £ 0 and q n = 1, 2, 3, ... 

For loads acting radially at the location r = T, 0, the general 
solution for the radial displacement, with q £ 0, q / -!p, is 

*< r > = 7 (1 - mfqT 1 ) (1 - r 2 ) - {1 + f sgn(r) 

- T r - | r | - tan(q) cos(q'r) cos(qr) 

(136) 

+ |sgn(r) sin(qr) cos(qr) - qta n(q) sin (^) sin (<l r ) 

+ i cos(qr) sin (q | r| ) , 

where 

-1, r < 0 

sgn(r) = 0, r = 0 (137) 

l, r > o . 


Equation 136 satisfies the boundary conditions and the transition condi- 
tions. Substitution of eq. 136 into eq. 125 yields the load-thrust 
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relationship. 


A£ 2 + Be + C = 0 


(138) 


where A, B, and C are functions of the thrust, X, and r". (See Appendix 
A). The solution procedure for the offset loading follows that of the 
midspan loading case. 


Eccentric Loading: q = 0 

For the thrust parameter equal to zero, the governing equation is 
again altered. The resulting displacement solution is. 


r 2 - r 

4 >(r) = 5 {4- (r + sgn (r ) ) - 4 (1 + r sgn(r)) 


+ ini t iu 2 

+ 12 " 


4 ~ sgn(r) + . 


(139) 


The loads for which the thrust is zero are given by the relationship. 




5 + 18 (f ) 2 + (r ) 4 
3-7 (r ) 2 + 5(T ) 4 -(r ) 6 


•]} - o , 


(140) 


which is obtained by substituting the functional form of the radial 
displacement into eq. 125. 


Eccentric Loading, q = — , n = 1, 2, 3, ... 

With q = n = 1, 2, 3, ... the adjoint problem outlined for the 
midspan loading must again be examined. The orthogonality condition 
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reveals that the applied loads for these values of thrust have the form, 

nit sin(-^) 

p= ^ . (141) 

cos(5fV 

For odd values of the radial displacement is. 


<t> n (r) = 7 (1 - r 2 ) + — 

q n 


_{1 + r sgn(r)} 

9 n Sin(q n ) cos(q n r) 


+ - ^ {7 + sgn(r)} + A* cos(qr) 

q n sin ( q n ) cos (q n r ) 


-g-i-— s 1n (q | r |) , 

V ln K) 


(142) 


while for even values of the displacement is given by. 


i ? i cos (q n r) * 

♦ n (r) -T (1 - r) +-4- (1 7 — ) * A n sin(q r) . (143) 

% cos (v 

★ 

The magnitude of A n for both displacement functions is found from the 
expression. 


A(A*) 2 + BA* + C = 0 , 


(144) 


where A, B, and C are constants which depend on the values of the 
thrust, and T (See Appendix A for the functional forms of A, B, and 

C). 
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Typical Numerical Results - Midspan Loading 

Shallow cylindrical ly curved panels exhibit a number of different 
deformation and buckling behaviors. The particular behavior depends on 
the parameter X, Some of the effects of orthotropic material properties 
on the value of X, and thus the behavior of curved panels, is explored 
in Appendix B. Generally, though, curved panels which have no eccen- 
tricity in the loading ff = 0) exhibit three distinct behaviors. In 
this section each type of behavior will be discussed in order to 
familiarize the reader with typical analytic results. These discussions 
will focus on the behavior under displacement control and load con- 
trol. However, the condition represented by the experimental setup to 
be described is displacement control. 

If X is sufficiently small the load-displacement behavior is mono- 
tonic and no snap-buckling occurs. In fig. 5 typical load-thrust, p vs. 
q, and load-midspan displacement, p vs. diagrams are shown for a 
panel which exhibits monotonic behavior. Here ^ denotes the 
nondimensional radial displacement at the location of the applied 
load. As the load and displacement increase from zero through points 1 
and 2, the stiffness (dp/dl>) of the panel decreases while the thrust 
increases. At point 3 the thrust is a maximum and the stiffness is a 
minimum. A further increase in load or displacement through points 4 
and 5, depending on which variable is controlled, causes the thrust to 
decrease and the stiffness to increase. Finally at point 6, a point 
close to the fully inverted configuration, the thrust is approaching 
zero. In fig. 6 the equilibrium configurations corresponding to the 
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FIG. 5 Load-thrust and load^ is placement 
RELATIONS, A =1.75, T =0 . 0 
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Note: Numbers correspond to points on fig, 5 



FIG, 6 Equilibrium configurations, \= 1 , 75 , r=0,0 
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FIG. 7 Surface strains vs. applied load, A =1.75, F =0.0 
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points denoted on the load-tnrust and load-displacement are shown. 

Figure 7 shows the behavior of the upper and lower surface strains, 
i.e., z = + -Tj- and z = - at the two quarter points (r = ± 0.5) and 
midspan (r = 0.0) , as a function of the applied load. Obviously due to 
symmetry of the response, the strains at the quarter points are 
identical. This panel exhibits no instability, a characteristic which 
is typical of panels with \ < 1.9762. 

A second type of behavior is known as limit point behavior. This 
is depicted in fig. 8. As the load, or displacement, is increased from 
zero through point 1, the thrust increases and the stiffness decrea- 
ses. At point 2 the load-thrust curve has a horizontal tangent and the 
stiffness of the panel is zero. At point 2 a relative maximum load, or 
limit point, is reached on the equilibrium path. Under load control 
this is a critical equilibrium configuration and increasing the load 
results in a sudden dynamic snap-through to a equilibrium configuration 
between points 5 and 6. Under displacement control this point is sta- 
ble. Increasing the displacement past point 2 will cause the load to 
decrease and the thrust to continue to increase. Point 3 denotes the 
maximum thrust for the panel. Increasing the displacement from point 3 
to point 4 causes both the thrust and the load to decrease. Notice that 
point 4 is also a load limit point. From point 4 through point 5 the 
load increases while the thrust decreases to zero with increasing dis- 
placement. Figure 9 shows the equilibrium configurations, under dis- 
placement control, for various values of load and thrust. Figure 10 
illustrates the strain behavior at the two quarter points and at mid- 
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FIG. 8 Load-thrust and loaded is placement 
RELATIONS, A =2.5, F=0.0 
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Mote: Numbers correspond to points on fig. 8 





A =2.5, r =0.0 


FIG. 9 Equilibrium configurations. 
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span. Under displacement control this panel exhibits no instability. 
Under load control there is the snap-through instability at point 2, the 
limit point, which is considered the buckling load. Limit point behav- 
ior is found for 1.9762 < \ < 2.825. 

The final type of behavior possible for the panel with no load 
eccentricity is known as bifurcation behavior. Bifurcation occurs on 
the initial loading path when the maximum thrust equals n. Figure 11 
shows typical load-thrust and load-displacement diagrams for bifurcation 
behavior. Load, displacement, and thrust all increase from no-load 
through point 1 to point 2. At point 2 a secondary equilibrium path, 
the straight line through points 2, 3, and 4, intersects the primary 
path. This secondary path represents the asymmetric response mentioned 
in Chapter 2 and included in eq. 127. (Actually there are two equilib- 
rium paths associated with asymmetric configurations intersecting at 
point 2 since two values of the amplitude b, eq. 130, are determined for 
each value of the load.) Point 2 is called a bifurcation point. In 
load control both the primary equilibrium path, beyond point 2, and the 
secondary equilibrium path configurations are unstable. Thus the panel 
will dynamically snap-through at the bifurcation point to a configura- 
tion on the primary path beyond point 6. Point 2 is considered the 
buckling load. In displacement control the secondary path is stable 
while the primary path becomes unstable. Increasing the displacement 
causes the load to decrease, along the straight line between points 2, 

3, and 4, and the thrust to remain constant. Note the load actually 
drops to zero at point 3. The spatial configuration of the panel 
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FIG, 11 Load-thrust and load-d is placement 
RELATIONS, A =4,0 , F =0.0 
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corresponding to points along the secondary path is asymmetric while the 
primary path corresponds to symmetric deformation. Point 4 represents 
tne second bifurcation point where the secondary path again intersects 
the primary path. Figure 12 illustrates the spatial deformation charac- 
teristics for a panel which exhibits bifurcation behavior and is under 
displacement control. Note the asymmetry in the deformed shape as the 
response moves along the bifurcated path. The surface strains shown in 
fig. 13 illustrate only those strains predicted for the stable equi- 
librium path under displacement control. 

Bifurcated solutions exist at all values of q = nit, n = 1, 2, 3, 

... n max , but only the first bifurcated path (n = 1) represents the 
physically important result for the static response. It should be noted 
here that bifurcation behavior is a mathematical result of an analysis 
for a perfectly symmetrical structure. Since, in practice, real struc- 
tures are not perfect this idealization will never be achieved. How- 
ever, with relatively small initial imperfections it could be closely 
approximated. Bifurcation behavior occurs for \ > 2.825. 

Typical Numerical Results - Eccentric Loading 

When offsets, or eccentricities, in the location of the applied 
loading are present, only two types of behavior exist. These are mono- 
tonic behavior and limit point behavior. Panels which exhibit monotonic 
or limit point behaviors when the load is at midspan show little change 
in their behavior for small eccentricities in the loading. Figures 14- 
16 show the monotonic response of a panel with geometric and material 
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Note; Numbers correspond to points on fig, 11 



FIG. 12 Equilibrium configurations, A =4.0, r = 0,0. 
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FIG. 13 Surface strains vs. applied load, A =/ 1 . 0 , r=0.0 
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Note: numbers correspond to points on fig. 14 
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FIG. 15 Equilibrium configurations, A =1.75, r =0-1 
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FIG. 16 Surface strains vs. applied load, A =1.75, T =0.1 
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properties identical to the panel of figs. 5-7, with the exception of a 
10% offset in the load being present. Here the load is offset to the 
right. The only significant effect of the offset is a slight asymmetry 
in the deformation. This is illustrated by the equilibrium configura- 
tions shown in fig. 15. Figure 16 shows the surface strain behavior for 
this loading case. 

Figures 17-19 show the effects of load offset on a panel which 
exhibits limit point behavior when the load is exactly at midspan. 

Again the only significant effect is a slight asymmetry in the deforma- 
tion. A comparison of figs. 8 and 17 shows that the offset leaves the 
limit load essentially unchanged. 

A panel which exhibits instability at bifurcation points when the 
load is exactly at midspan is affected to a much greater degree by 
eccentricities in the load position. This is easily seen by comparing 
the load-thrust and load-displacement diagrams of fig. 11, the case of 
midspan loading, with those of fig. 20. Following the equilibrium path 
shown in fig. 20, an increase in the displacement causes the thrust and 
load to increase up to point 2. Here point 2 represents a load limit 
point (horizontal tangent). For the midspan loading case, fig. 11, 
point 2 was not a limit point. It was a bifurcation point. For the 
offset loading, any increase in the applied load beyond point 2 would 
cause loss of stability. On the other hand, controlling the displace- 
ment from point 2 to point 4 causes the load to decrease and the thrust 
to continue increasing. At point 4, q > it. Thus for load offsets q 
= i is no longer the limiting value of thrust for the stable equilibrium 
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FIG, 17 Load-thrust and load-displacement 
RELATIONS, \ =2.5 , T =0.1 
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Note: numbers correspond to points on fig. 17 






FIG. 13 Equilibrium configurations, \ =2.5, r = 0.1 
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FIG. 19 Surface strains vs. applied load, X =2 ' 5 ^r=0-1 
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path. Recall from fig. 11, q = n was the limiting value for the midspan 
case. This allows the eccentrically loaded panel to exhibit higher 
strain levels than are found with midspan loading. Also, at point 4 a 
displacement limit point (vertical tangent) is found. Therefore an 
increase in the displacement causes the arch to snap-through to point 
5. The load level p suddenly jumps, as does the thrust level. From 
point 5 an increase in the displacement causes the load to increase 
while the thrust decreases to zero. A comparison of figs. 20 and 11 
reveals that the bifurcation load of fig. 11 is significantly higher 
than the limit load of the offset response in fig. 20. Figure 21 shows 
the equilibrium configurations for various points along the equilibrium 
path. Of special interest are the shapes associated with points 4 and 
5. The shape corresponding to point 4 represents the equilibrium con- 
figuration just before the displacement limit point while the shape 
corresponding to point 5 is the configuration of the panel just after 
the displacement limit point. Figure 22 illustrates the behavior of the 
surface strains vs. the applied load. It should be noted here that the 
sharp jumps in the strain curves are due to the displacement limit 
point. This portion of the curve is not stable, the surface strains 
experience a dynamic change due to the change in load and thrust. 

From the analysis of the response, it appears that the behavior of 
a cylindrical panel is quite sensitive to load eccentricities for cer- 
tain ranges of X. To fully understand the effect of load offsets, the 
relation between the magnitude of the buckling load and the offset angle 
is examined. Since the problem is physically symmetric with respect to 
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Note: Numbers correspond to points on fig. 20 



FIG. 21 Equilibrium configurations, A = ^-0, P=0.1 
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FIG. 23 Sensitivity of buckling load to load eccentricity 
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the midspan vertical, only positive load offsets are shown. Figure 23a 
shows the effect of load eccentricity on the buckling load for three 
different panels, all of which exhibit limit point behavior when the 
load is exactly at midspan. Here the limit load is considered the 
buckling load. For panel 1, which has k = 1.99, increasing the offset 
angle causes only slight increases in the buckling load. However, with 
this panel, an offset causes the behavior to switch from limit point to 
monotonic behavior. Thus there no longer is an instability. Panel 2, 
with k = 2.05, resembles panel 1 with the exception that limit point 
behavior applies for the entire range of load offsets. Panel 3, 
where k = 2.5, is different than the first two in that the buckling load 
initially decreases with increasing load offset. Still the panel is not 
overly sensitive to eccentricities in the loading. Small perturba- 
tions (T « 0.1) in the load position only cause negligible changes in 
the buckling load. Both panels 2 and 3 exhibit buckling loads approach- 
ing infinity as the applied load is moved very close to the supports. 

A panel which exhibits bifurcation instability when the load is 
exactly at midspan tends to be quite sensitive to load eccentricities. 
Figure 23b shows the buckling load as a function of offset angle for 
such a panel. The bifurcation load corresponds to the buckling load 
for 7 = 0, and the limit point load corresponds to the buckling load for 
all other load offsets. The buckling load shows significant reductions 
even for small load offsets. This type of behavior leads to the panel 
being considered imperfection-sensitive. 


Chapter 4 

EXPERIMENTAL SET-UP AND PROCEDURE 

Little experimental work has been performed in the area of the 
large, out-of-plane deformation behavior of shallow composite cylin- 
ders. The work that has been presented tends to be somewhat limited in 
scope in that the range of geometry is small or the effects of load 
eccentricities have not been investigated. To check the validity of the 
analysis just presented, a series of tests on composite specimens of 
varying geometries and load offsets were performed. These tests were to 
determine if all the behavioral characteristics predicted by the analy- 
sis are observable. Also, the effectiveness of the solution procedure, 
as well as classical lamination theory, in determining the large defor- 
mation response of composite structures was to be verified. 

In order to make a comparison of the analytic and experimental 
results, the conditions of the experiment needed to closely approximate 
those used for the analysis. This required the test specimen to be 
supported such that the ends remained free of rotational constraints 
while the span between supports was kept constant. Also the load needed 
to act radially at the panel middle surface. To meet these require- 
ments, as well as others, a test fixture and test procedure were 
designed to facilitate a fair comparison of theory and experiment. 

Curved Panel Specimens 

The test specimens were fabricated by NASA-Langley using AS4/3502 
graphite-epoxy pre-preg. Table 1 shows the specified geometry and 
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Table 1 


Intended Test Specimen Geometry 


Specimen 

Number 

Lay up 

Thickness 

(in.) 

Load 

Offset 

(%) 

1 

[(90/0) 4 ] s 

0.08 

0 

2 

C(90/0) 4 ] s 

0.08 

15 

3 

[(90/0) 4 ] s 

0.08 

30 

4 

C(90/0) 4 ] s 

0.08 

50 

5 

[(90/0) 3 ] s 

0.06 

0 

6 

C(90/0) 3 ] s 

0.06 

15 

7 

[ ( 90/0 ) 3 ] s 

0.06 

30 

8 

E(90/0) 3 ] s 

0.06 

50 

9 

[90/45/0/-45] s 

0.04 

15 

10 

[90/45/0/-45] s 

0.04 

30 

11 

[90/45/0/-45] s 

0.04 

50 

12 

[(90/0) 4 ] s 

0.08 

0 

13 

E(90/0) 4 ] s 

0.08 

15 

14 

C(90/0) 4 ] s 

0.08 

30 

15 

C(90/0)] s 

0.02 

0 

16 

C(90/0)] s 

0.02 

15 

17 

C(90/0)] s 

0.02 

30 

18 

C(90/0)] s 

0.02 

50 

Note: For 

all cases the intended radius 

was 60 in. with 

an arc length 

of 

12 in., with the exception of 

specimens 12-14, 

where the arc 


length was 6 in. 
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stacking sequence for each specimen tested. For each geometry which was 
to be tested, a large 18 x 24 in. curved panel was fabricated using 
existing curved forms. After curing, each large panel was ultrasonic- 
ally C-scanned to determine if any irregularities were present. For 
each large panel the results of the C-scan showed uniformity in the 
transmittance of sound waves, suggesting a uniformity of the material. 
Upon determining the quality of the panels, the specimens were then 
machined from the larger panels. Each specimen had a set of three holes 
drilled at the point of load application. This will be discussed 
shortly. By keeping a tight tolerance on the arc-wise location of the 
holes, it was possible to control the amount of load eccentricity. 

After machining, slotted dowels 0.375 in. in diameter, to be further 
discussed in connection with the test fixture, were attached to each end 
of the specimens. Figure 24 shows 6 in. arc length and 12 in. arc 
length panels both with and without the dowels attached. Two dowels are 
shown in fig. 24, as are two larger dowels which will be discussed. 
Strain gaging has been attached to the panels in fig. 24 and a 6 in. 
rule is visible in the photograph. 

Test Fixture 

Figure 25 shows a schematic of the test fixture while fig. 26 shows 
a photograph of the fixture. The major components of the fixture 
are: a) end dowels, b) support blocks, c) a base plate, d) loading 

head and pivot, and e) tongue and connector. Each component and its 
function will now be described. 
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FIG. 24 Photograph of machined specimens 






FIG. 25 Schematic of test fixture 
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FIG. 26 Photograph of test fixture 
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a) End Dowels 

The end dowels were 0.375 in. diameter steel rods, slotted to 
accommodate the ends of the curved panels. The dowels were one part of 
a system designed to accurately represent a simple support condition. 

The ends of the panel were potted solidly in the dowel slots using 
epoxy. The epoxy prevented the ends of the panel from brooming under 
the inplane compressive thrust loads. This potting is shown in fig. 

27. The rods were slotted along a rod radius and the panels were cen- 
tered, using shims if necessary, within the slots. With care, the end 
of the middle surface of the panel could be made to coincide with the 
centerline of the dowel. With this arrangement, the boundary conditions 
at the ends of the middle surface (S = ±S 0 ) could be accurately 
enforced. 

b ) Support Blocks 

The support blocks were an important component in the system and 
two designs were used before the test fixture was satisfactory. The 
support blocks provided the simple support condition by preventing 
translational displacements but allowing rotational displacement. The 
original design called for support blocks into which semicircular 
grooves were machined. The grooves were designed to accommodate the 
circular dowels. A lubricant between the dowel and the groove was used 
to minimize frictional effects related to rotation of the ends of the 
specimen. Unfortunately, due to the inplane thrust forces induced in 
the curved panel as it deformed, significant frictional forces developed 
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END DOWEL 


PANEL 



FIG. 27 Schematic of end dowels 
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in the end blocks. These frictional forces resulted in end moments on 
the panel, a violation of the moment-free boundary condition assumed in 
the analysis. A second design of the end block used six roller bearings 
in each block. These roller bearings were situated so that the dowels 
were forced against them. This is shown in fig. 28. There were 
actually three bearings in a line above the end dowel and three in a 
line below the end dowel. Figure 29 shows a photograph of the roller 
bearing end blocks. 

c ) Base Plate 

Figure 30 shows a schematic of the base plate. The base plate is 
also visible in figs. 25 and 26. The function of the base plate was to 
provide stiffness to the entire fixture and to provide a solid founda- 
tion for locating the support blocks. The base plate was 19 x 4 x 1 in. 
and was machined from a larger piece of steel stock. Care was taken to 
insure that the top and bottom surfaces of the plate were flat and 
parallel with each other. One support block was mounted at the location 
of the six isolated holes (top of fig. 30) while the other was mounted 
in either set of slots. Hardened steel pins, 0.250 in. in diameter and 
protruding from the base plate into the support block, were also used to 
secure the block mounted at the holes. More will be said of this 
shortly. The two sets of slots were to accommodate the panels of two 
different lengths. The slots were machined into the plate so that the 
position of the support blocks relative to one another could be 
adjusted. This was necessary since the panels all had slightly differ- 


FIG. 28 Schematic of support block 


ROLLER BEARINGS 


END DOWEL 


ROLLER BEARINGS 


HOLE FOR ATTACHMENT 
TO BASE PLATE 



FIG. 29 Photograph of support block 
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FIG, 30 Schematic of base plate 
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ent lengths when machined. The block at the slotted end could be 
adjusted along the length of the base plate so that the panel was tight 
in the fixture but without the radius of curvature being too different 
than the natural value, i.e., the value specified in Table 1. This was 
very important since the response of curved panels can be quite 
sensitive to the initial radius of curvature. 

There was concern that the large inplane forces generated as the 
panels deformed would cause the blocks to move relative to the base 
plate. This is the reason hardened steel pins were used with the one 
support block. For the support block attached at the slots, another 
design was initially incorporated. A third steel block was bolted into 
position at the six holes near the slots. This block was designed to be 
a back-up support for the block at the slots. The load was transferred 
to the back-up by way of several horizontally oriented bolts jammed 
between the back-up block and the support block. As it turned out, 
there was no slipping between the base plate and the support blocks at 
either end of the panel and so the back-up block was not used. 

d) Loading Head and Pivot 

The analysis assumed that the panel was subjected to a line load 
applied at the geometric middle surface. There were two difficulties in 
achieving this ideal. First, it was physically impossible to contact 
the geometric middle surface. A line load acting radially from above 
would contact the outer surface of the panel. This would result in a 
slight offset between the location the load was assumed to act and where 
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it actually did. This is not serious until the cross section of the 
panel at the line load location begins to rotate. Then this effect, 
coupled with the load, would induce a small moment in the panel, in 
addition to the line load. This could have been accounted for in the 
analysis but such a moment was deformation-dependent, making the govern- 
ing equations more difficult to solve. Thus to have the line load act 
at the geometric middle surface, a novel loading mechanism was 
adapted. Incorporated into this mechanism was a scheme to insure 
against the second problem, namely that no moment be transferred into 
the panel through the loading fixture. Figure 31 shows the loading head 
and pivot fork used to transmit a pure line load to the geometric middle 
surface of the panel. Figure 32 shows the pivoting action of the mecha- 
nism. Three small bolts through the panel and the large dowels (shown 
in fig. 24) secured the arc-wise location of the line load. Centering 
the panel in the gap between the large dowels guaranteed that the line 
of pivoting passed through the geometric middle surface. Of course the 
pivots themselves prevented any moment from being transferred to the 
panel. Figure 33 shows a photograph of the pivot. 

e) Tongue and Connector 

The pivot fork was attached to the bottom end of two long steel 
plates. The long plate assembly, or tongue, is shown in fig. 34. In 
the analysis it was assumed the line load moved downward radially and 
deflected the panel, with the simply supported ends of the panel remain- 
ing stationary. In the experiment, the vertical displacement of the 
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THROUGH GAP 


FIG. 31 Schematic of loading head and pivot 


PIVOT FORK 




FIG. 32 Illustration of pivoting action 
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FIG. 33 


Photograph of loading head and pivot 
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CONNECTOR 


TONGUE 



PIVOT FORK 


FIG. 34 Schematic of tongue and connector 
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panel at the location of the line load was held to zero and the simply 
supported ends moved upward. Due to the assummed shallowness of the 
panels tested this was essentially the same as holding the radial 
displacement to zero at this point. The tongue was used to attach the 
panel to the stationary upper head of the load frame. The base plate 
rested on the lower movable head. To deform the panel, the lower head 
was moved upward. The tongue was made as long as possible for the 
following reason: As the panel deformed, tangential (u) displacements 

resulted. This meant that the line load position moved tangentially. 

In the analysis it was assumed that the line load acted strictly in a 
radial direction. To do this the line load must translate tangentially, 
the line of action remaining parallel to the original no-load 
direction. Experimentally it was very difficult to translate the line 
load tangentially. An infinitely long tongue, fixed at the top and 
attached to the panel at the bottom, would satisfy this requirement. 
Practically speaking, a long tongue is all that was possible. The 
maximum tongue length was dictated by the physical dimensions of the 
load frame. A connector simply attached the upper end of the tongue to 
a universal joint and ultimately the fixed upper head. Figure 35 shows 
another overall view of the test fixture with a panel in place. 

Instrumentation 

After fabrication the specimens were instrumented with 350 ohm 
strain gages with a gage length of 0.25 in. The excitation voltage for 
the strain gages was 2 volts. The strain gages were mounted back-to- 
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back on the specimens so that the middle surface strain could be moni- 
tored. The location of the gages was determined by the amount of load 
offset for the specimen. For load offsets of 15% of the semi-opening 
angle and less, the gages were placed at the quarter points. When the 
load offset exceeded 15% the gages were placed at midspan and at the 
quarter point furthest from the load. To measure the radial displace- 
ment, dial gages mounted radially were employed. The dial gages were 
capable of measuring displacements of 0.001 in. and greater. Again the 
location of the dial gages depended on the degree of the load offset. 

For offsets of 15% and less, the dial gages were located at the quarter 
points and at the point of load application. For offsets greater than 
15% the dial gages were placed at the point of load application, mid- 
span, and at the quarter point farthest from the load. The tests were 
carried out in the Engineering Science and Mechanics department of 
Virginia Tech using a displacement-controlled Instron loading frame. 

The load was determined using a 1000 pound Instron load cell placed in 
series with the load fixture. The strain gage data was acquired using a 
system consisting of Vishay 2120 Wheatstone bridges and signal condi- 
tioner, an HP3495A scanner with an A/D converter, and a Tektronix 4051 
computer. The dial gage readings were manually entered into the Tektro- 
nix computer while strain gage and load data were automatically recorded 
on digital tape. Figures 36 and 37 shows the experimental set-up with 
instrumentation in place. 



FIG, 36 Photograph of experimental set-up 

side VIEW 
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FIG. 37 Photograph of experimental set-up 


END VIEW 
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Initial Measurements 

Before the actual testing of the specimens took place, some initial 
measurements related to panel geometry and material properties were 
taken. The arc length, width, and thickness were measured to assure 
each was as specified. Also the positioning of the holes drilled for 
load positioning were checked for accuracy. Tensile tests to determine 
the extensional stiffness, and three point bend tests to determine the 
bending stiffness were made. The tensile tests were performed on strips 
cut from the larger panels from which the specimens were cut. The bend 
tests were conducted on the actual specimens. Details of the test 
procedure and results are contained in Appendix C. 

The initial radius of each specimen was measured after it was 
fitted into the test fixture. Measuring the radius after the panel was 
in the fixture assured that the measured radius was indeed the initial 
radius of the specimen tested. To accomplish the measurement of the 
radius, a dial gage and a traversing mechanism, capable of measuring arc 
length location to within 0.05 in., were used to obtain the spatial 
coordinates of several points on the mounted specimen. By using a 
parabolic least-square fit to the data, the average radius of the 
specimen was determined. Table 2 shows the measured dimensions for the 
specimens tested. 

Test Procedure 

Having made the initial measurements on the specimen, excluding the 
radius measurement, it was placed loosely in the test fixture. In this 
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Table 2 


Measured Test Specimen Geometry 


Specimen 

Number 

Radius 

(in.) 

Thickness 

(in.) 

1 

56.34 

.084 

2 

52.04 

.082 

3 

54.72 

.083 

4 

51.00 

.082 

5 

55.42 

.060 

6 

55.34 

.062 

7 

50.35 

.063 

8 

50.09 

.063 

9 

52.56 

.042 

10 

55.41 

.043 

11 

57.33 

.042 

12 

64.7 

.082 

13 

62.4 

.081 

14 

63.8 

.081 

15 

50.41 

.022 

16 

— 

-- 

17 

-- 

— 

18 




Note: Of specimens 15-18 only specimen 15 was actually tested 

thus specimens 16-18 had no initial measurements made. 
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position the strain gages were connected and they were balanced and 

zeroed. The specimen was then tightened in the fixture by sliding the 

movable support block. To assure a tight fit the specimen was slightly 
compressed by the support blocks until a slight bending strain was 
registered. The loading head was then attached and the levelness of the 
fixture and the vertical orientation of the tongue were checked and 
corrected as needed. The initial radius of curvature measurements were 
then taken. After this, the load cell reading was zeroed. The dis- 
placement rate was set at 0.01 in. /min. and readings of load, strain, 
and displacement were recorded at every 0.01 in. of midspan displace- 
ment. To record the dial gage readings, the loading head displacement 

was stopped. The readings were then entered into the Tektronix. During 
the first tests the support blocks were monitored with dial gages for 
any sign of elastic or rigid body displacement. Since the blocks showed 
no sign of bending or slipping relative to the base plate, this was 
discontinued for later tests. 


Chapter 5 

CORRELATION OF EXPERIMENTAL AND THEORETICAL RESULTS 

In order to validate an analysis it is necessary to compare it with 
experimental results. Here the analytic results presented in Chapter 3 
are compared to experimental results obtained using the apparatus and 
test fixture described in Chapter 4. So as to make the comparison with 
the experiment clearer, the analytic results shown represent only those 
configurations assumed to be stable. The initial measurements made on 
the test specimens concerning the geometry and material properties were 
used in generating the analytical curves presented with experimental 
results. Also, the small initial strain in the panels, due to fitting 
the specimens into the fixture, has been subtracted from the experi- 
mental data presented. In all tests performed this initial strain was 
less than 50pe in magnitude, or less than 2.5% of the maximum strain. 
Both the experimental and analytic results are presented in nondimen- 
sional form, with the exception of the shape data, where the units are 
inches. 

In general, the agreement between theoretical and experimental 
results was good. All the trends predicted by the analysis were 
observed in the experimental results. Although the agreement was seen 
to be good, there was still some deviation between theory and experi- 
ment. This deviation may be attributed to a number of factors. Perhaps 
the most significant was the fact that the specimens were found to have 
imperfect shapes, i.e., they were not exactly circular. The work done 
by Schreyer [4] shows that geometric imperfections have much the same 
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effect on panel response as do load eccentricities. This is the most 
likely reason for not seeing a distinct bifurcation of the equilibrium 
path when the load was exactly at midspan. Other possible factors for 
deviation of the experimental results from the predictions include the 
possibility of friction-generated moments at the supports, even though 
great care was taken to minimize this. The friction was manifested by 
small load jumps, characteristic of end support sticking, observed 
during the testing of the panels. Thickness variations observed in the 
test specimens during the initial measurements were also likely to cause 
differences between experiment and theory. This thickness variation 
could be responsible for the differences between the predicted and 
observed maximum surface strains seen for some of the test specimens. 

In this section the experimental results from selected test cases 
are presented. These cases, along with their associated parameters, are 
given in Table 3. Since the trends observed for most of the panels 
tested were similar, the test cases shown here are representati ve. 
Results from the [90/0] s panels, which were quite thin, are not shown 
here. Problems with these specimens developed during the testing. It 
appeared that friction generated at the supports was, at times, 
sufficient to constrain free rotation of the specimen. The specimen 
would actually bend at the support ends and so the specimen behaved at 
times as if it had clamped ends instead of simply supported ends. 

Figures 38-40 show the predicted and observed response for a 
[(90/0)4]$ com P osite specimen with a 12 in. arc length, a 56.34 in. 
radius, and no load offset, test case 1. The analytic predictions are 
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Table 3 

Test Cases Presented 


Test 

Case 

Number 

Specimen 

Number 

Lay-up 

Radius 

(in.) 

Arc 

Length 

(in.) 

Offset 

Angle 

(%) 

1 

1 

C(90/0) 4 ] s 

56.34 

12 

0 

2 

2 

to 

1 — 1 

o 

o 

CT> 

1 — 1 

52.04 

12 

15 

3 

3 

to 

1 — 1 

O 

0 
(T> 

1 1 

54.72 

12 

30 

4 

5 

C(90/0] 3 ] s 

55.42 

12 

0 

5 

6 

1 — 1 

to 

o 

0 

GO 

1 1 

(/> 

55.34 

12 

15 

6 

9 

[90/ 45/0/ -45] s 

52.56 

12 

15 

7 

12 

C(90/0) 4 ] s 

64.7 

6 

0 


0 . 021 ” 
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AND LOAD-DISPLACEMENT RELATIONS; R=56 14 , 8=12.0,, 

r-O.O, I <90/0) 4 1 s 
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FIG. 39 Comparison of theory and experiment: equilibrium 
configurations; R=56,34, S=12 . 0 , T =0 . 0, 
C(9Q/0) 4 1 s 
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r = -0.5 




FIG. 40 Comparison of theory and experiment: surface 
STRAINS VS. APPLIED LOAD; R=56,34, S=12,0, 

r=0.0, C(90/0 ) 4 ] s 
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represented by the solid line while the open symbols represent the 
experimental data. The analysis predicts that such a panel will exhibit 
bifurcation of the primary equilibrium branch. The load-displacement 
diagram of fig. 38 shows that the observed behavior closely follows the 
predicted response. (The darkened symbols in fig. 38 correspond to 
decreasing values of thrust.) From the discussion of bifurcation 
behavior in Chapter 3 this suggests that the initial imperfections 
associated with the panel are relatively small. The load-thrust diagram 
of fig. 38 shows more clearly that the actual behavior is limit point 
behavior. Since it has been shown that geometric imperfections effect 
the panel response in much the same way as do load eccentricities, it 
would be beneficial to refer to the discussion in Chapter 3 of the 
effects of load eccentricities on panels exhibiting bifurcation behavior 
when they are loaded exactly at midspan. Appendix D investigates the 
effect of geometric imperfections on the panel behavior. It should be 
noted here that the inplane stress resultant, N 0 , is assumed to be 
constant (see eq. 80). Perhaps, this is not exactly true and the devia- 
tion of theory and experiment is due to a slight spatial variation of 
the thrust. To experimentally obtain the thrust parameter, q, the back- 
to-back strain gages at each spatial position are averaged to obtain the 
middle surface strain. The thrust parameter is then calculated, using 
eqs. 65, 105, and 123, at each point and is then averaged. Despite the 
symmetry to the set-up, for a given specimen and load, the thrust para- 
meter at each strain gage location did not have the same numerical 
value. This suggests that q may not have been spatially constant. 
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Alternatively, a slight misalignment of the strain gages or other errors 
associated with strain gage measurements could have made a very accurate 
determination of q impossible. Also, it is seen from the load-displace- 
ment diagram that the initial panel stiffness, represented by the slope 
of the equilibrium path, observed experimentally differed from that 
which was predicted. This difference could be attributed to errors made 
in the material property and geometry measurements of the specimens. 

The parameter x, which governs the panel response, is sensitive to these 
measurements. Appendix E explores the effect of the initial panel 
stiffness on the overall response. Fig. 39 shows the panels geometric 
configuration of various displacement levels. The panel configurations 
predicted by the analysis are purely symmetric up to the bifurcation 
point. After this the configurations become asymmetric. Though not 
shown in fig. 39, during the experiment it was observed that the 
configurations began to show asymmetry at much lower load levels than 
were predicted. This again led one to believe that geometric 
imperfections were present. Figure 40 shows a comparison between 
theoretical and measured strains. The maximum surface strains predicted 
by the theory are somewhat higher than the observed strains. This may 
have been due to local thickness variations, since only an average panel 
thickness was used in the analysis. 

Figures 41-43 show the response for a [(90/0)4] s 12 in. arc length 
composite specimen with a 52.04 in. radius and a 15% load offset, test 
case 2. With no load offset such a panel shows bifurcation response. 
Ideally, it would have been desirable to have a panel with a radius of 
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FIG. 42 Comparison of theory and experiment: equilibrium 
CONFIGURATIONS; R=52,04, S=12 . 0, F =0 . 15, 

C (90/0) 4 ] s 


102 




FIG. 43 Comparison of theory and experiment: surface 

STRAINS VS. APPLIED LOAD; R=52.04, S=12 . 0, r = 0 • 15, 
C(90/0) 4 3 s 
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curvature equal to the first experimental case described. Since the 
initial radius of curvature was quite sensitive to the span between 
support blocks, it was difficult to obtain the same radius for a group 
of specimens. However, since radii of curvature for the first and 
second experimental cases were quite close, and since the two radii of 
curvature would produce bifurcation behavior in the no-offset case, it 
is fair to say that the only difference between the first and second 
experimental cases was the 15% load offset. This second experimental 
case performed, then, as expected. A comparison of this offset load 
case and the previous midspan load case shows a decrease in the experi- 
mentally determined limit loads, thus verifying the prediction that such 
a structure is sensitive to load eccentricity. The load-thrust and 
load-displacement diagrams of fig. 41 show good agreement between theory 
and experiment. The same is true of the shape predictions. The surface 
strain predictions, shown in fig. 43, deviate from the observed strains 
significantly. 

Figures 44-46 show the results for a [(90/0)4] $ panel with a 30% 
offset. The panel arc length was 12 in. and the radius of curvature was 
54.72 in. The comparison between theory and experiment for the load- 
thrust relation was not as good as for the other two first cases but 
still the basic trends were there. For offset loadings, the maximum 
thrust tends to increase with increasing load offset for panels which 
exhibit bifurcation behavior when the load is exactly at midspan. Thus 
the force on the supports increases. During the testing of this speci- 
men it was observed that friction at the supports had more of an effect 
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FIG, 45 Comparison of theory and experiment: Equilibrium 
configurations; R=54.72, S=12.0, P =0.3, 
[(90/0) 4 ] $ 
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FIG. 46 Comparison of theory and experiment: surface 

STRAINS VS. APPLIED LOAD; R=54.72, S=12.0, F =0,3, 
C(90/0) 4 ] $ 
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than for the previous tests. The friction was detected from slight load 
jumps due to intermittent sticking of the supports, thus constraining 
rotation. The configurations predicted, fig. 45, show fair agreement 
with the shapes observed. The predicted and observed strains once again 
deviate from one another. 

Figures 47-49 show the response for a [(90/0)3] $ 12 in. arc length 
panel with a radius of 55.42 in. and no load offset while figs. 50-52 
show the response of a very similar panel (55.34 in. radius) with a 15% 
load offset. For the midspan loading the correlation between theory and 
experiment is good, particularly for the load-displacement relation and 
the strains. Again there appeared to be some geometric imperfections 
present since limit point behavior, rather than bifurcation behavior, 
was observed for the specimen with no load offset. This was also born 
out by the observation of the equilibrium configurations during test- 
ing. The shapes observed were asymmetric in nature at much lower load 
levels than were predicted. In general, however, shape correlation was 
acceptable, as can be seen in fig. 48. The offset loading case also 
shows extraordinary agreement between theory and experiment. Load- 
deflection, load-thrust, shape, and the surface strain predictions 
correlate well with the observed results. 

The results for a [90/45/0/-45] s 12 in. arc length panel, with a 
radius of 52.56 in. and a 15% load offset, test case 6, are shown in 
fig. 53-55. In order to theoretically predict the response of a quasi- 
isotropic composite, it was assumed that the response was one-dimen- 
sional and that the only nonzero curvature change was in the 0 direc- 
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FIG. 47 Comparison of theory and experiment: load-thrust 
AND LOAD-DISPLACEMENT RELATIONS; R=55,42, S=12,0, 

r=0.0, [(90/0) 3] s 
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FIG. 48 Comparison of theory and experiment: Equilibrium 
configurations; R=55.42, S=12.0, F =0,0, 

C(90/0) 3 1 S 
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FIG. 49 Comparison of theory and experiment: surface 

STRAINS VS. APPLIED LOAD; R=55,42., S=12 . 0, F =0 . 0 , 
C ( 90/0 ) 3 ] s 
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FIG. 50 Comparison of theory and experiment:, load-thrust 
AND LOAD-DISPLACEMENT RELATIONS; R=55.34, S=12,0, 

r-o.is, [00/0)3] s 






FIG. 51 Comparison of theory and experiment: Equilibrium 
configurations; R=55,3 4, S=12.0^ r = 0.15^ 

[(90/0) 3 ] s 
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FIG. 52’ Comparison of theory and experiment: surface 
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AND LOAD-DISPLACEMENT RELATIONS; R=52.55, $=12.0, 
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FIG. 54 Comparison of theory and experiment - : equilibrium 
configurations; R=52.56, S=12 . 0 , r-0.15, 
[90/45/0/-45] s 
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FIG, 55 Comparison of theory and experiment: surface 
STRAINS VS. APPLIED LOAD; R=52,56, S=12.0, 

r=0.15, C 90/45/07-^5] s 
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tion. Thus a reduced bending stiffness was used. In reality, D 16 and 
D 25 being nonzero invalidate this assumption. Due to the constraint 
placed on the test specimen by the fixture it was thought that these 
one-dimensional assumptions were a good first approximation. A compari- 
son of theory and experiment shows good agreement. The correlation 
between the observed and predicted load-displacement, deformed shapes, 
and strains are also seen to be quite good. 

Figures 56-57 show the response of a [(90/0)^] s 6 in. arc length 
panel with a radius of 64.7 and no load offset. This is case 7 and is 
the last case discussed here. The predicted response is limit point 
behavior but the observed response is that of monotonic behavior. The 
reason for this discrepancy between theory and experiment is not 
immediately obvious. Flexibility of the supports is not an issue here 
since the forces exerted on the support blocks are less than in previous 
cases. Error in the initial radius measurement has also been ruled out 
since it was double-checked and found to be correct. Also, the initial 
material property measurements, recalling their effect on the 
parameter \, are felt to be accurate enough as to not be a factor 
here. It seems likely that there is a "start-up" problem with the 
experimental set-up. By "start-up" is meant low load and thrust levels, 
i.e., initial portion of test. This could be the result of mechanical 
slack or dead zones in the fixture. Recall that for the previous test 
cases there was a discrepancy between the predicted and observed initial 
stiffnesses. Appendix E ruled out the possibility of error in initial 
measurements as the cause of this discrepancy. Since the initial mid- 
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FIG, 56 Comparison of theory and experiment: load-thrust 
AND LOAD-DISPLACEMENT RELATIONS; R=64,7, S=6,0, 

r-0.0, C(90/0) 4 1 s 
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FIG. 57 Comparison of theory and experiment: surface 

STRAINS VS, APPLIED LOAD; R=64,7, S=6,0,F=0.0, 
C(90/0) 4 3 s 
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span rise of the panel of test case 7 is significantly smaller than for 
the other test cases, it is possible that the start up problem effects a 
much greater portion of panel response, thus creating the poor correla- 
tion between prediction and observation. Whatever the cause, it was 
observed in tests of the same panel with load offsets of 15% and 30%. 
Thus the phenomenon was repeatable. 

In order to address the question of load eccentricity, the limit 
point loads for the various load offsets for a given group of specimens, 
i.e., same lay-up and arc length, were plotted along with an upper and 
lower bound on the theoretical predictions. Here the upper bound corre- 
sponds to the curve for the panel with the smallest value of \, while 
the lower bound corresponds to the curve for the panel with the largest 
value of \. Since the radius of each panel within the group changes, it 
is necessary to plot the upper and lower bound. Figure 58 shows the 
plot of limit load vs. percent load offset for the [(90/0)^] s 12 in. arc 
length specimens. From this figure it seems that the sensitivity to 
load eccentricities predicted is indeed close to that which was 
observed. Similarly fig. 59 shows the sensitivity plot for the 
C ( 90/0 ) 3 ] s 12 in. arc length specimens, and fig. 60 shows the sensi- 
tivity for the [90/45/0/-45] s specimens. It appears the observed sensi- 
tivity has been accurately predicted. 


CD 


121 



. 58 Comparison of theory and experiment: sensitivity 
of buckling load to load eccentricity, [(90/0) 4 ]g 
specimen group 
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FIG. 59 Comparison of theory and experiment: sensitivity 

OF BUCKLING LOAD TO LOAD ECCENTRICITY, C(90/0) 3] $ 
SPECIMEN GROUP 
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FIG. 60 Comparison of theory and experiment: sensitivity 

OF BUCKLING LOAD TO LOAD ECCENTRICITY, 

[90/45/0/-45] $ specimen group 


Chapter 6 

CONCLUSIONS AND RECOMMENDATIONS 


The study presented in this thesis was intended to extend the 
knowledge of the static behavior of shallow orthotropic curved panels. 
The principle aims of the study were to determine the effects of ortho- 
tropic material properties on the behavior of shallow curved panels, to 
experimentally verify the analysis, and thus ultimately lay the ground 
work for a dynamic analysis of the same problem. An exact solution, 
taking into account load eccentricity, was developed for a simply 
supported circularly cylindrical curved panel subjected to a line load 
along a cylinder generator. Also, an experimental apparatus was 
designed to load the curved panels. The analytical predictions were 
compared with experimental measurements. 

From this study the following conclusions may be made: 

1. Analytically, three types of panel behavior were seen for 
panels loaded exactly at midspan; monotonic, limit point, and 
bifurcation. Only two types of behavior were exhibited when 
the load was offset; monotonic and limit point. The type of 
behavior predicted depends solely on one parameter, \. Geome- 
try, material properties, and stacking sequence influence \ for 
an orthotropic panel. For an isotropic panel \ is affected by 
geometry only. 

2. Those panels which exhibit bifurcation behavior when loaded 
exactly at midspan were seen to be imperfection-sensitive. 

This was demonstrated both analytically and experimentally. 
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3. Initial geometric imperfections were detected and quantified 
while making initial radius measurements. Due to the imperfec- 
tions present, inherent of any real structure, bifurcation 
behavior was not seen. However the observed response did 
closely approximate bifurcation behavior. 

4. The overall correlation between theory and experiment was 
good. This implies that the solution method used is reasonably 
accurate. Predictions of the limit loads were quite accurate. 

5. The predictions of strain in the panels were fair for the most 
part. However, a spanwise spatial variation in the middle 
surface strain was observed during testing. This is counter to 
the assumption that the strain was constant along the span. 

6. Due to the small bending stiffness inherent to thin laminates, 
it was impossible to obtain accurate results for these 
laminates with the current experimental set-up. The friction- 
related moments generated at the ends of the simply supported 
panels at times caused these panels to behave as if they were 
clamped. 

7. In analyzing the test results there seems to have been "start- 
up" or initial load and deformation problems with the experi- 
mental set-up. This is borne out by the lack of correlation 
between prediction and observation for the low thrust 
regimes. This would explain the discrepancies in initial 
stiffnesses as well as the total lack of correlation for test 
case 7 in which the maximum thrust was expected to be rela- 

ti vely small . 
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As a result of this work, a number of topics can be singled out for 
further investigations. Of these topics, a number are relatively minor 
extensions of the previous work. These are: 

1. Modify the analysis to include different boundary conditions; 
e.g., frictional effects. 

2. Include the effects of geometric imperfections in the analysis. 

3. Conduct further static tests with various material properties 
and geometries. By varying the material properties the effects 
of orthotropicity described in this text may be experimentally 
verified. 

4. Investigate the "start-up" problem encountered during testing. 

In addition to the above, recommendations requiring significant 
extensions of this work are: 

1. Extend the analysis to include the effects of shear deformation 
and to take into account angle ply laminates. This would allow 
a far greater number of laminates to be investigated and would 
perhaps be an intermediate step to the dynamic analysis. 

2. Extend the analysis into the dynamic regime. Ultimately it is 
hoped that an understanding of the behavior of a composite 
fuselage, upon impacting a flat, rigid surface, may be 
obtained. By doing so the crashworthiness of a composite 
fuselage may be assessed. 
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Appendix A 

COEFFICIENTS OF LOAD-THRUST RELATIONSHIPS 


The solution technique applied to the governing differential equa- 
tions discussed in Chapter 3 requires the solution of a quadratic equa- 
tion. This equation relates either the applied load or a constant of 
the general solution of <j> (r ) to the thrust. Here the functional forms 
of the coefficients to the quadratic equation are presented. For each 
case the coefficients, found to be transcendental functions of the 
thrust, are found by substituting the functional form of the radial 
displacement into the integrated strain-displacement relationship, 
expressed nondimensional ly by eq. 125. 

Midspan Loading: q ^ 0 and q f n = 1, 3, 5, ... 

The load-thrust relationship is of the form, 


AC 2 + Be + C = 0 


(Al) 


when q f n-rt, n = 1, 2, 3 


Here A, B, and C are. 



4 4 5 

16q^ 16q^ 16q 3 


(A2) 



(A3) 



_ i + i si n(2q) 


(A4) 
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When q = nit, n = 1, 2, 3, ... the integrated strain-displacement rela- 
tionship yields, 

A£ 2 + B£ + C = Db 2 . (A5) 

This relationship relates the applied load to the magnitude of the 
asymmetric displacement. The coefficients are. 


A ■ - 3 4 

16q^ 


(A6) 

B ■ [1 

- (-1)"] 

(A7) 


1 

(A8) 

ii 

Q 


(A9) 
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Midspan Loading: q = — , n = 1, 3, 5, ... 

The integrated strain-displacement relationship yields. 


A(A*) 2 + BA* t C = 0 , 


(A10) 


to determine a constant, A n> from the general solution of <t> n (r) given by 
eq. 134. Here the coefficients are. 


A = 


4q‘ 


(All) 
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4 ( - i)n 

(A12) 

- q2 1 + 3 + 1 / n n 

7 '< ‘ 

( A 1 3 ) 


Eccentric Loading: q ^ 0 and q ^ — , n = 1, 2, 3, ... 

The load-thrust relationship, when load eccentricity is present is, 

AC 2 + Be + C = 0 . (A14) 


For the offset load case the coefficients of the quadratic equation are 
functions of the load offset as well as the thrust. The coefficients 
are: 


3 _ (r ) 2 _ r sin (2qr ) _ sin 2 (2qr ) 

16q 4 8q 4 16q 4 sin(2q) 16q^tan(2q) 


2 _ 

, cos (qr )sin(2q) 

+ 5 

16q 3 

o 

3cos (qr )tan(q) 
16q 5 


sin * 2 (qr ) 


4q xan(2q)sin(2q) 


sin 2 (2qr )tan(q) sin 4 (qr ) 


16q' 


8q°tan (q) 


(A15) 


+ 5sin 2 (qr) 
16q 5 tan(q) 


B _ r*sin (qr) tanqcos (qr) + cos (qr) 1 

4q 3 cos(q) 4q 3 cos(q) 2q 4 cos(q) 2q 4 

r = iL _ I _ sin(2g ) + 1 

X 4 6 8q 3 cos 2 (q) 4q 2 cos 2 (q) 


(A 16 ) 
(A17) 
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Eccentric Loading: q = y-» n = 1, 2, 3, 

For both odd and even values of n, the integrated strain-displace- 
ment relationship yields. 


A(A*) 2 + BA*+ C = 0. 
x n ' n 


(A18) 


For odd values of n the coefficients are: 




(A 19 ) 


B = £°. s . 2 (q r ) (_i) n ( f _i) . sin 2 (gr ) ( _ 1} n 

q q 


+ 7fq < _1 ) n E cos (2qr ) - COS ( 2q ) ] 


(A20) 
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For even values of n the coefficients are: 


A = - - 3 - 

M 4 


(A22) 


B = 0 


(A23) 


C .1.51 + _2_ 

6 X 4 4q 2 


(A24) 


Appendix B 

EFFECT OF ORTHOTROPY ON PANEL BEHAVIOR 
The parameter \ is defined as. 


where k22> The extensional stiffness, and D£ 2 > The bending stiffness, 
are given by. 


Concentrating on cross-ply laminates, where each lamina is composed of 
the same material, the reduced transformed stiffness is. 


n 



(B2) 


(B3) 


Q 22 » f° r 0° laminae 



(B4) 


Q^, for 90° laminae. 


The Q's are the reduced stiffness and are defined as. 



(B5) 



(B6) 
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where is the Young's modulus in the fiber direction and E 22 is the 
Young's modulus perpendicular to the fiber direction. By definition it 
is assumed that > E 22 * 

Using the above definitions for the stiffnesses, a general expres- 
sion for the extensional stiffness and the bending stiffness can be 
developed in terms of the number of laminae, the stacking sequences, and 
the elastic moduli of the laminae. For convenience these expressions 
are shown only for the [(90/0) n ] s family of laminates. For these lami- 
nates, the extensional and bending stiffnesses are, 


A 


22 


r < E n- e ?2> 

1-v 12 v 21 


(B7) 


D 


22 " 


N 2 t 3 


[-(N+3)E n + (N-3)E 22 ] 
1 - v 12 v 21 


(B8) 


where N is the number of laminae (N=4n) and t is the thickness of each 
lamina. Substituting eqs. B7 and B8 into eq. B1 yields. 


" = 4 N V[(1 ♦ 

n L 22 


* )] r2 


(89) 


where use is made of the fact S = 2Rp. 

To emphasize the importance of taking the orthotropy of the mater- 
ial into account, three shallow circularly cylindrical panels are com- 
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pared as to the behavior exhibited. Each panel is assumed to have the 
same geometry, only the material used in the fabrication of the panels 
Is varied. For the comparison aluminum, glass-epoxy, and graphite-epoxy 
panels are used. Table B1 summarizes the geometry, material properties, 
and calculated values of X for each of the panels. The value of X for 
the aluminum panel is such that bifurcation on the primary equilibrium 
path would occur when the panel is loaded exactly at midspan. However, 
when glass-epoxy or graphite-epoxy is used, limit point behavior would 
be expected. Thus by increasing the degree of orthotropy, the panel 
behavior changes from bifurcation, where a asymetric response would be 
expected when in displacement control, to limit point behavior, where a 
jump in load and thrust would occur. 

For a composite material, the number of laminae, as well as the 
degree of orthotropy, effect the value of X, Figure B1 shows the effect 
of the number of laminae, as well as the degree of orthotropy present, 
on the calculated value of X, With an increasing number of laminae, the 
effect of orthotropy becomes less pronounced. 

Stacking arrangement, though not explored specifically here, could 
be important in determining x and thus panel response. Though they 
would have the same weight, two panels each with 16 laminae could be 
made to behave quite differently, depending on how the laminae are 
arranged. 
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Table B1 

Panel Geometries and Material Properties 

Panel 



Aluminum 

G1 ass-Epoxy 

Graphite-Epoxy 

Arc Length (in.) 

6 

6 

6 

Radius (in.) 

45 

45 

45 

Thickness (in.) 

0.08 

0.08 

0.08 

Lay up 

— 

[90 4 /0 4 ] s 

[90 4 /0 4 ] s 

Ell (P si ) 

10 x 10 6 

7.8 x 10 6 

25 x 10 6 

E 22 (psi) 

10 x 10 6 

2.6 x 10 6 

1 x 10 6 

e h/e 22 

1.0 

3.0 

25.0 

X 

2.94 

2.72 

2.58 



K90/0) n 3 $ WITH CONSTANT 
LAMINATE THICKNESS 
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FIG. B1 Effect of orthotropy and lay-up on panel response 
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Appendix C 


DETERMINATION OF EXTENSIONAL STIFFNESS AND 
BENDING STIFFNESS 

As mentioned in Chapter 4, the initial measurements made on the 
test specimens included the determination of the extensional and bending 
stiffnesses. To accurately compare the results of the analysis and the 
experiment it is necessary to have correct material properties to input 
into the analysis. Here the test set-ups and procedures used for the 
measurement of both the extensional stiffness and bending stiffness are 
descri bed. 

Extensional Stiffness 

The extensional stiffnesses for each layup were obtained from 
tensile test data. A tensile test specimen was machined from each of 
the large panels, as were the curved panel test specimens. The tensile 
coupon, 0.5 x 12 in., was cut from the panel such that the stiffness in 
the x direction was actually measured. Since the specimens were cross- 
ply laminates, the stiffnesses in the x direction and e direction should 
have been equal. The specimens were tested in the Instron Tensile 
Machine using an Instron load cell to measure applied load. Back-to- 
back strain gages were mounted at the center of the gage section. The 
strain gages used were 120 ohm gages with approximately 4 feet of lead 
wire. The excitation voltage for the strain gages was 2 volts. The 
strain gage and load data were acquired by using a digital data acquisi- 
tion system. For each test performed the slope of the load vs. strain 
curve was calculated using a least squares curve fit. This slope. 
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proportional to the laminate Young's modulus in the x direction, was 
then divided by the width of the tensile specimen to obtain the 
extensional stiffness. Each test was performed twice in order to obtain 
results as accurately as possible. The results for the layups investi- 
gated are summarized in Table Cl. 

Bending Stiffness 

The bending stiffnesses for each of the curved panel specimens 
tested were obtained by performing three-point bend tests on the speci- 
mens themselves. The specimens were placed on the fixture so that the 
panel was simply supported, restricting end displacement in the vertical 
direction while allowing freedom to displace horizontally. The load was 
applied through a hanger located at midspan to which calibrated dead 
weights were attached. A schematic of the test fixture used for the 
three-point bend tests is shown in fig. Cl. A dial gage capable of 
resolving displacements to 0.001 in. was used to measure the vertical 
displacement at midspan. Load vs. displacement data was recorded and 
the slope determined using a least-squares fit to the data. Assuming a 
linear straight beam theory was suitable for describing the behavior of 
the bend specimen, the bending stiffness was obtained from the 
expression. 


D = 

U 22 48w 


(Cl) 
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Table Cl 

Measured Panel Stiffnesses 

Test No. Lay-up A 22 (lb/in) D 22 (lb-in) 


1 

CO 

1 1 

o 

o 

1 — 1 

7.8194 x 10 5 

616.2 

2 

[(90/0) 4 ] s 

7.8194 x 10 5 

606.3 

3 

C(90/0) 4 ] s 

7.8194 x 10 5 

617.2 

4 

C(90/0) 4 ] s 

7.8194 x 10 5 

612.7 

5 

C(90/0) 3 ] s 

5.85 x 10 5 

263.5 

6 

C(90/0) 3 ] s 

5.85 x 10 5 

265.4 

7 

C (90/ 0 ) 3 ] s 

5.85 x 10 5 

268.2 

8 

C(90/0) 3 ] s 

5.85 x 10 5 

260.1 

9 

[90/45/0/-45] s 

3.175 x 10 5 

90.9 

10 

[90/45/0/-45] s 

3.175 x 10 5 

90.5 

11 

[90/45/0/-45] s 

3.175 x 10 5 

93.3 

12 

C(90/0) 4 ] s 

7.8194 x 10 5 

613.5 

13 

C(90/0) 4 ] s 

7.8194 x 10 5 

609.8 

14 

C(90/0) 4 ] s 

7.8194 x 10 5 

614.2 

15 

[90/0] $ 

2.21 x 10 5 

12.36 

16 

1 1 

LO 

o 

o 
1 — 1 

CO 

2.21 x 10 5 

12.70 

17 

l — l 
iO 

o 

0 

1 1 

CO 

2.21 x 10 5 

12.30 

18 

[90/0] s 

2.21 x 10 5 

12.51 


DIAL GAGE 



FIG. Cl Three-point bend test fixture 
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where p is the load intensity, w is the midspan displacement and L is 
the span between supports. The results of the three-point bend tests 
are summarized in Table Cl. 


Appendix D 

EFFECT OF GEOMETRIC IMPERFECTIONS ON PANEL BEHAVIOR 

As was mentioned earlier, Schreyer [4] found the effects of geo- 
metric imperfections on the behavior of isotropic curved panels to be 
qualitatively the same as the effects of load offsets. The geometric 
imperfections investigated in Schreyer' s study were of the same form as 
the asymmetric deformation mode of the curved panels. In the present 
study it was thought that similar geometric imperfections were partially 
responsible for some of the deviation between theoretical and experi- 
mental results given in Chapter 5. Since the location of the holes 
drilled into the panel for locating the loading head were measured and 
felt to be accurately centered, load eccentricity was not felt to be a 
problem. 

Figures 38-40 show the predicted and observed response for a 
[ (90/0)43 s composite specimen with no load offset. The behavior expec- 
ted of this panel was bifurcation behavior. While the load-displacement 
diagram of fig. 38 appears to show bifurcation behavior, the load-thrust 
diagram leads one to believe that some imperfection may be present. For 
this test case a number of measurements were made along the entire span 
of the specimen for determining the average radius. (For subsequent 
tests these measurements were made only over half the span.) From these 
measurements it was possible to determine the deviation of the initial 
configuration from the assumed perfect cylindrical configuration. 

Figure D1 shows the measured deviation of the initial configuration 
from a perfect cylindrical geometry. The perfect configuration, y 0 (x). 
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FIG. D1 Measured initial geometric imperfection 
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is assumed to be a cylinder with a radius of curvature, -p which is 
determined by a least-squares fit of the data taken along the span of 
the curved panel. The measured configuration is denoted by y (x ) • From 
this figure it is seen that the deviation, or imperfection, tends to 
take on the form of the asymmetric deformation mode discussed earlier. 
Since geometric imperfections were found from the initial measurements 
and were found to be of the same form as the asymmetric deformation 
mode, as was the case in Schreyer's work, it was a distinct possibility 
that some of the differences between theory and experiment are due to 
these imperfections. 

To explore this point further a 5% load offset was compared with 
the theoretical midspan loading and the experimental points. The 5% 
load offset was used to show the qualitative effect of a small imperfec- 
tion in the initial configuration. Figure D2 shows the theoretical 
curves, both a 0% and a 5% load offset, and experimental data for the 
load-thrust and load-displacement relations for the panel of test case 1 
(see Chapter 5). The 5% offset case tends to more closely approximate 
the data taken from the test over a greater portion of the response 
curves. 

Finally, fig. D3 shows the spatial shapes for the first 18 load or 
displacement points of fig. D2. The predicted shapes are of course 
symmetric with respect to midspan. For the first 6 or 7 data points, 
the measured shapes are also symmetric. However, the 8th point shows 
that the actual shape has lost symmetry, the left side deflecting more 
than the right side. At points 13 and 14 the asymmetry has shifted, the 
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relations; R=56.34, S=12.0, r=0.0 and 
F=0.05, [(90/0) 4 ] s 
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FIG. D3 Comparison of theory and experiment: Equilibrium 
CONF IGURAT IONS j R=56.34, $=12 . 0, F =0 . 0, 

[(90,0) 4 ] $ 
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FIG. D3 (cont.) 
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FIG. D3 (cont.) 
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right side deflecting more. These deviations in symmetry of the spatial 
shape are felt to be due directly to initial imperfection of fig. Dl. 

Even though the work presented in this appendix is only for one 
specific case, in all likelihood some form of geometric imperfection is 
present in each of the other test cases as well. The results shown here 
show that the geometric imperfection found was of the appropriate form 
to cause an imperfect response and that by taking the initial imperfect 
configuration into account, the predicted response more closely models 
the panel behavior. 


Appendix E 


EFFECT OF INITIAL STIFFNESS ON PANEL BEHAVIOR 


In Chapter 5 it was noted that the initial stiffness, i.e., the 
initial slope of the load-displacement relation, predicted by the theory 
and the observed initial stiffness did not correlate well. Perhaps some 
of the discrepancy between the theory and experiment was due to this 
initial slope difference. Here the effect of the initial stiffness on 
the overall response is explored with the expectation that, by correctly 
modeling the linear behavior, the nonlinear behavior can be closely 
modeled. The initial panel stiffness may be determined by examining the 
linear response of the panel. In order to determine this stiffness, the 
theorems of Castigliano are used for the no-load offset case. A full 
investigation of this is found in [11]. 

The complimentary strain energy for the panel is. 


For a shallow panel the statically admissible forces, found from a free 
body diagram, are 


where T is the horizontal support reaction. Substituting these forces 
into eq. El yields 



(El) 



(E2) 



(E3) 
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22 


(E4) 


in which 


. °22 2R 2 p 4 


(E5) 


22 15 



(E6) 



(E7) 


Applying Casti gl iano 's second theorem to the complimentary strain energy 
determines the horizontal reaction in terms of the applied load. Then, 
by applying Castigliano's first theorem to the complimentary strain 
energy an expression for the initial stiffness is obtained. The initial 
stiffness is. 


Rp 4C T c p - C tp 


(E8) 


Substituting for the constants Cj, Cyp, and Cp the stiffness is 
rewritten as. 


K 



(E9) 


1 375 ( 

7 " TTS [ 



where use was made of the relationships. 


S = 2R(3 


(E10) 
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2„4 


a 22 r p 

~ D 


22 


(Ell) 


From eq. E9 it is seen that the geometric and material properties 
of the panel govern the initial response. In order to model the non- 
linear response accurately it seems that the predicted linear response 
should match the observed response. This may be accomplished by 'jugg- 
ling' the panel parameters in the theoretical formulation until the 
predicted response matches that of the observed. Table El shows the 
effect of changing the panel parameters on the initial stiffness of the 
panel. Each column shows the effect of a percentage decrease of the 
parameters. The initial parameters used are those of test case 1 from 
Chapter 5. 

Upon examining the predicted and observed initial stiffnesses it is 
seen that the parameters must be changed substantially to match the 
stiffness. By changing the panel parameters so dramatically, the pre- 
dicted character of the nonlinear response is changed considerably. 

This is shown in fig. El. In this figure, the initial stiffness in the 
analysis is adjusted, through \, to match the observed initial 
stiffness. It is seen that although the initial portions of the load- 
displacement and load-thrust agree, the overall comparison is very 
poor. Thus it seems inaccuracies in initial stiffness could not be 
responsible for the differences between experiment and theory observed 
in this study. Also, since small inaccuracies, say 5%, in the panels 
geometric and material properties have such a small impact on the 


154 


Table El 

Effect of Panel Parameters on Initial Stiffness 


K/K 0 



1% 

5% 

10% 

R 

1.008 

1.037 

1.076 

°22 

.994 

.968 

.935 

A 22 

.997 

.981 

.961 

S 

1.015 

1.077 

1.152 


K o = K ( R o* D 22 q » A 22 q » S o) 
R 0 = 56.34 in. 

Do? = 616.2 lb-in. 
cc o 

A?? = 7.8194 x 10 5 lb/in. 
cc o 

S Q = 12.0 in. 
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FIG. El Comparison of theory and experiment (with 

- 0.02 MATCHED .INITIAL STIFFNESS) : R=56,3 i !, 

S=12.0, r=0.0, t(90/0) ] s 
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initial stiffness, it seems unlikely that poor measurements of the 
panel's geometric and material properties would be responsible for 
differences between observed and predicted panel response. 
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